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A Proposal for improving the Hadronization of QCD currents in 
TAUOLA 

P. Roig a * 

a Instituto de Fisica Corpuscular, IFIC, CSIC-Universitat de Valencia. 
Apt. de Correus 22085, E-46071 Valencia, Spain 

After overviewing the general features of semileptonic decays of the tau lepton, I will recall the most widely used 
model for them, namely that of Kiihn-Santamaria (KS), and I will explain the subsequent works that were done 
along these lines and that are implemented in the TAUOLA library for analysing tau decays. After that, I will 
move to the description of our project that aims to achieve a theory as close as possible to QCD for the considered 
decays. I conclude by emphasizing the importance of the implementation of our work in TAUOLA. 
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1. INTRODUCTION 

Hadronic decays of the r lepton are neatly split- 
ted into the electroweak part and the strong one, 
responsible for the hadronization. While the first 
one is well under control, the second one still 
eludes our understanding from first principles pQ . 
The division of both sectors yields a clean envi- 
ronment for studying the hadronization proper- 
ties of QCD currents. 

The decay amplitude for the considered pro- 
cesses can be written as: 
Gf 

M = --j= VcKMUv T l»{l - l^UrH^ , (1) 

where the hadronic vector, 7i M , encodes our lack 
of knowledge of the precise hadronization mecha- 
nism: 

= ( {P( Pl )}ti I (v„ - z lCQCD lo > • (2) 

Symmetries help us to decompose 7i M depending 
on the number of final-state mesons, n. In the 
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three meson case this reads: 

W M = V llM Ff(Q 2 ,s u s 2 ) + V 2/M F^(Q 2 , s u s 2 ) + 

Q^(Q 2 , Sl ,s 2 )+ iV 3 ^(Q 2 , Sl ,s 2 ) , (3) 



and 



V 2fl = [g^v - 



(P2-PI) V , 



~W 

= (Pi + P2 + P3)fi ,Si = (Q- Pif . (4) 

Here F{, i = 1,2,3, correspond to the axial- 
vector current {A^) while drives the vector 
current (V^). The form factors F\ and F 2 have 
a transverse structure in the total hadron mo- 
menta, Q' 1 , and drive a J p = 1 + transition. The 
scalar form factor, i 7 ^, vanishes with the mass of 
the Goldstone bosons (chiral limit) and, accord- 
ingly, gives a tiny contribution. 

It is clear that one cannot derive the F{ from 
QCD. The aim of this communication is to con- 
vince that the most adequate way out is the use 
of a phenomenologically motivated theory that, 
in the energy region spanned by tau decays, re- 
sembles QCD as much as possible. First works 
on the considered modes were [2] and [3]. 

In any case, the approximate symmetries of 
QCD are useful. They rule what is the theory to 
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be used in its very low-energy domain and guide 
the construction of higher-energy theories, as I 
will sketch in section [31 Before of that I will ex- 
plain the model implemented in the MonteCarlo 
generators that is presently used in these decays. 

2. KS MODEL IN TAUOLA 

The Kiihn-Santamaria model [2] proposes phc- 
nomcnologically motivated amplitudes that, by 
construction, fulfill the LO {0{p 2 )) xPT result 
in the chiral limit. However, it fails to reproduce 
the NLO result [J], [5J- It relies on the success- 
ful vector meson dominance [B] and ensures the 
right behaviour at high energies, guaranteed by 
the Breit-Wigner propagators accounting for res- 
onance exchange. Though widely used, its link 
with QCD is not proved and the proposed off-shell 
widths are not QFT-motivated. In any case, it 
provides a good description for the two and three 
pion decays of the tau. 

There is a lot of literature following this work 
studying many other hadronic decays of the r 
0, [E], M- On the other hand, TAUOLA has 
grown over the years jTD] to be a complete li- 
brary providing final state with full topology in- 
cluding neutrinos, resonances and lighter mesons 
and complete spin structure throughout the de- 
cay. In these works, the hadronization part of 
the matrix elements follows the KS model. In the 
remaining of the section I will focus on [7] describ- 
ing the decays r — > KKkv t , t — > Kttkv t . In it, 
the chiral, one- and two-resonance exchange con- 
tributions are obtained by expanding the follow- 
ing products of Breit-Wigner's with correspond- 
ing factors weighting the relative contribution of 
each resonance: 

M ? Mj 

Mf - s - i^Ti (s) Ml - t - iV~tT 2 (i) ' 

Throughout our work [11] , we have noticed 
that these studies lack from several contributions 
that happen to be non-negligible in our formal- 
ism, namely, they only account for p 0, ~ and K*° 
being exchanged in one channel in the modes 
K+K'ir-, K-K°ir°, K-Tr-ir + and xVtt-. 
That is, they are only including these exchanges 
either in V/ 1 or ijl}, while one naturally ob- 



tains them in both. The authors of Refs. [7] and 
[5] needed to include noticeable different masses 
and widths for the p' resonance in the vector and 
axial-vector currents. Another point we have no- 
ticed is the inclusion of a different number of mul- 
tiplets in the vector and axial-vector form factors 
(3 and 2, respectively). In summary, we may con- 
clude that although the KS model describes quite 
well the two and three pion decay modes of the 
t, it seems that its generalization for other three 
meson modes needs improvement 0. The target of 
our work, whose guidelines I will be explaining in 
the following, is to tackle an ambitious program 
for providing a description as close as possible to 
QCD for the considered decays. 

3. QCD IN TERMS OF MESONS 

The key idea is to pursue a description in terms 
of the relevant degrees of freedom in the range of 
energies spanned by tau decays that incorporates 
the approximate symmetries of QCD therein. 
Well below the p(770) mass, xPT [12], [13] is the 
effective theory of QCD. Nevertheless, it only ap- 
plies to a small part of the phase space in r decays 
[T4] . so one needs to include the resonances as 
explicit degrees of freedom. One expects those of 
spin one to be the most important ones, according 
to vector meson dominance 6 . This emerges as 
a result in Resonance Chiral Theory (RxT) [15] . 
which is built upon the approximate chiral sym- 
metry of low-energy QCD for the lightest pseu- 
doscalar mesons and unitary symmetry for the 
resonances being the general procedure guaran- 
teed by Weinberg's Theorem [12], [16]. In this 
case, writing a phenomcnological Lagrangian in 
terms of mesons fulfilling the basic symmetries 
of light-flavoured QCD will provide a right de- 
scription. The perturbative expansion of RxT is 
guided by its large- N c limit [T7]. At LO in 1/N C , 
meson dynamics is described by tree level dia- 
grams obtained from an effective local Lagrangian 
including the interactions among an infinite num- 
ber of stable resonances. We include finite reso- 
nance widths (a NLO effect in this expansion) 
within our framework 18J. We are also depart- 
ing from the strict Nc — > oo limit because we 

2 Current experimental situation is reviewed in section [5] 
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consider just one multiplet of resonances per set 
of quantum numbers (single resonance approxi- 
mation |19| ) and not the infinite tower predicted 
that cannot be included in a model independent 
way. For convenience, we choose to represent the 
spin-one mesons in the antisymmetric tensor for- 
malism |20) . 

The relevant part of the RxT Lagrangian is 

m- m\- m, m- 

C RxT = + X+ > + ^(Wf > 

( V^vTu" ) + {^{A, v r ) + ci n 

5 7 

+£kin + 2-s ^*®VAP + T7~" OvJP 
t=l <=1 v 

4 5 

+ £<fcO», VP + £^0< WJ ,, (6) 

i— 1 i—l 

where all couplings are real, being .F the pion de- 
cay constant in the chiral limit. The notation 
is that of Ref. [TS]. P stands for the lightest 
pseudoscalar mesons and A and V for the (axial)- 
vector mesons. Furthermore, all couplings in the 
second line are defined to be dimensionless. For 
the explicit form of the operators in the last line, 
see [21], [22], [23]. 

In order to inherit the maximum possible fea- 
tures of QCD, to be implemented in RxT, we 
still have to exploit the matching between or- 
der parameters of spontaneous chiral symmetry 
breaking with partonic QCD related quantities. 
The matching of n-point Green Functions in the 
OPE of QCD and in RxT has been shown to be 
a fruitful procedure [5D], [H], [22], [21], E3, [H]. 
Additionally, we will demand to the vector and 
axial-vector form factors a Brodsky-Lepage-like 
behaviour [27] , 

While symmetry fully determines the struc- 
ture of the operators, it is the QCD-ruled short- 
distance behaviour who restricts certain combi- 
nations of couplings rendering RxT predictive 
provided we fix these few remaining parameters 
restoring to phenomenology: Fy could be ex- 
tracted from the measured r(p° — > e + e~), Gy 
from r(p° — ► 7r + 7r~), Fa from r(ai —> irj) and 
the Ai's from T(a± — > prr) which stars on the 



t — » 3ttv t processes themselves. T(uj — > 777), 
T(u> — * 3tt) and the 0(p 6 ) correction to T(tt — ► 
77) may give us information on the remaining 
couplings [22J. However, we are not going to use 
these values but to constrain the couplings by de- 
manding the right short-distance behaviour to the 
form factors appearing in r decays as we will see 
in section [4] 

4. ASYMPTOTIC BEHAVIOUR AND 
QCD CONSTRAINTS 

There are 24 unknown couplings in £r x t © that 
may appear in the calculation of three meson de- 
cays of the t. They got reduced when computing 
the Fcynman diagrams involved. In addition to 
Fy and Gy, there only appear three combina- 
tions of the {Aj}| =1 , four of the {ci}J =1 , two of 
the {dj}| =1 and four of the {<7i}f =1 - The number 
of free parameters has been reduced from 24 to 
15. 

We require the form factors of the and 
currents into KKtt modes vanish at infinite trans- 
fer of momentum. As a result, we obtain con- 
straints |28] among all axial-vector current cou- 
plings but Ao, that are also the most general ones 
satisfying the demanded asymptotic behaviour in 
r — ► 3niy T . Proceeding analogously with the vec- 
tor current form factor results in five additional 
restrictions [28]. From the 24 initially free cou- 
plings in Eq. ([6]), only five remain free: C4, 
ci + c 2 + 8c 3 - c 5 , di + 8d 2 - d 3 , 54 and 55. 
After fitting T(u> — » 3rr) -using some of the rela- 
tions in Refs. [22], [29] - only C4 and (74 remain 
unknown. 

Relating the measured isovector component 
of e + e" — * KKtt to the total isovector cross- 
section to KKtt 30 and employing CVC to 
relate the latter to r — > KKtt [31] we were 
able 32J to provide a theoretical expression for 
<t (e + e~ — > KKtt) that we fitted to BaBar data 
[33] obtaining c 4 = -0.052 ± 0.003 and 54 = 
— 0.20jZo'ii S Remarkably, this procedure has 
allowed us to be sensitive to the sign of C4 im- 

3 However, we realized later that the relation used in |30) 
to relate the measured isovector component in e+e~ — > 
KgK^ir^ 1 to e+e — — > KKn is not correct. For a detailed 
explanation on this issue, see 
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proving the work done in [28 . Our results for 
the decay widths of the considered channels are 
consistent with the PDG values [Mj. We plan to 
confront our predictions for the spectra with the 
forthcoming experimental data. 

5. EXPERIMENTAL SITUATION 

CLEO [35], BaBar [36] and Belle [37] have been 
collecting good quality data onr-* KKirv T de- 
cays. CLEO [35 announced that the parameter- 
ization in Ref. [7] was unable to fit their data. 
They modified it by including two additional pa- 
rameters in order to obtain a good fit, but at the 
prize of violating a property of the strong inter- 
actions, namely the Wess-Zumino [38] normaliza- 
tion emanating from the chiral anomaly of QCD, 
as was put forward in Ref. [3S] • The conclusion to 
be drawn from their work is clear: CLEO stated 
that the KS model could not fit the data and the 
alternative parameterization they proposed is for- 
bidden from first principles. We would like that 
our expressions were used to analyse forthcom- 
ing data on the three meson decays of the r. Al- 
though BaBar has not published these results yet, 
they have managed [33j to split with great pre- 
cission the isoscalar and isovector contributions 
in e + e~ — > KKtt so that, under CVC we can use 
it for the considered r decays, as it has been done 
in Ref. [31]. We expect eagerly both BaBar and 
Belle new data to settle the issue of an adequate 
description of these decays. 

6. CONCLUSIONS 

I have reviewed the main assumptions of the KS 
model and its generalizations. We have found in- 
consistencies in the latter where the former did 
not fall into. Namely, the inclusion of a particle 
(p 1 ) with different properties depending on where 
it appears, a different spectra of propagating reso- 
nances (2 or 3 multiplets) in the vector and axial- 
vector current form factors or the lack of several 
non-negligible contributions in some channels. 

I have explained the phenomcnological ap- 
proach that we follow that includes every rele- 
vant (and known) piece of QCD: We preserve the 
correct chiral limit at low energies that gives the 



right normalization to our form factors, we em- 
ploy large- Nq QCD arguments to use our theory 
-predicting vector meson dominance- in terms of 
mesons and it fulfills the high-energy conditions 
of the fundamental theory at the mesonic level. 
While we remain predictive by including only the 
lightest multiplet of axial-vector and vector reso- 
nances, our Lagrangian can easily be extended in 
a systematic way to account for the contribution 
of higher resonance states, whose parameters may 
be fitted to experiment. 

Theory and experiment are linked through 
event generators. The precision achieved in theo- 
retical computations and the accuracy sought in 
the measurements justifies the effort that is be- 
ing made in improving the MonteCarlo's. Con- 
cerning the hadronic matrix elements appearing 
in r decays, an estimate of the error has tradi- 
tionally been obtained by comparing the results 
yielded by implementing the Brcit-Wigner fac- 
tors used in the KS model ([5]) with those in- 
corporating resonance exchange a la Gounaris- 
Sakurai [40] , The inclusion of our expressions 
for the hadronic matrix elements in the TAUOLA 
[TO] library for r decays and in PHOKHARA gT] 
for the hadronic e + e~ cross-section will be useful 
in order to describe better the hadronization of 
QCD currents. In the aforementioned Working 
Group there was unanimity in that understand- 
ing the hadronic contribution in both kinds of 
processes is mandatory for the MonteCarlo's em- 
ployed at the r-charm- and B-factories. This task 
will be of great help also in LHC and any future 
collider. 
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